arXiv:gr-qc/0211038vl 12 Nov 2002 


COSMIC EVOLUTION IN 
GENERALISED BRANS-DICKE THEORY 


B.K.Sahoo 1 and L.P.Singh 2 
Department of Physics, 

Utkal University,Bhubaneswar-751004,India 
1 bijaya@iopb.res.in 
2 lambodar@iopb.res.in 


Abstract 

We have studied the Generalised Brans-Dicke theory and obtained 
exact solutions of a(t ),cj>(t), and cu(t) for different epochs of the cosmic 
evolution. We discuss how inflation,decceleration,cosmic acceleration 
can result from this solution.The time variation of G(t) is also exam¬ 
ined. 
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1 INTRODUCTION 

The Brans-Dicke Theory(BD) is defined by a scalar field 0 and a constant 
coupling function a;[1].This is perhaps the most natural extension of gen¬ 
eral theory of relativity (GR) which is obtained in the limit of uo —» oo and 
0= constant[2],This theory interestingly appears naturally in supergravity 
theory,Kaluza-Klein theories and in all the known effective string actions. One 
important property of BD theory is that it gives simple expanding solutions[ 
3]for field 0(f) and scale factor a(t) which are compatible with solar sys¬ 
tem experiments[4],The Generalised Brans-Dicke theory(GBD)[2],sometimes 
referred to as graviton-dilaton or scalar -tensor(ST) theory, is instead, de¬ 
fined by tv which is an arbitrary function of the scalar field 0 (dilaton) and 
is thus implicitly a function of time ca(0(f)).Hence,it includes a number of 
models,for every function u. Naturally,a few attempts have been taken to 
study the dynamics of the universe using this formalism.Some of the cos¬ 
mological problems studied under this formalism includes nonsingular uni¬ 
verse^],amplification of gravitational wave[6],solutions for different epochs of 
cosmic evolution and time variation of Newton’s gravitational constant G[7] 
etc.These problems are studied for a specific forms of u;(0) involving quite 
complicated time-dependence for a(t),0(f).Thus,it is difficult to draw a clear 
physical picture about the evolution of the universe. 

In order to carry out a detailed study of the dynamics of the cosmic evo¬ 
lution in this formalism,knowledge about exact time-dependence of a(t),0(f) 
and u{t) for different epochs is essential. In a previous work[8],we have 
obtained power law time-dependence of cn(f) with the power of time deter¬ 
mined in terms of the exponent of scale factor a(t) which is taken to vary 
as a(f) = ao (j~) -To fix the value a we had taken the help of observational 
information.Similar result is also arrived at by Diaz-Rivera and Pimentel 
[9].However,in these works the value of a can not be fixed purely by way of 
solving the gravitational and scalar field equations. In the present work we 
take our previous attempt a little further and explore to see if a consistent 
solution of the GBD equations of motion based on simple power-law tempo¬ 
ral behaviour of a,0 and c o can fix the time exponents. 

We have arranged this paper in the following manner.In the section 2 we 
obtain the solutions for a,0,cn for different epochs.We present a discussion 
about the present value of u in section 3 and examined the time variation of 
Newton’s gravitational constant G in section 4 and conclude with section 5. 
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2 FIELD EQUATIONS AND SOLUTIONS 

For a universe filled with perfect fluid and described by Friedman- Robertson- 
Walker space-time with scale factor a(t) with spatial curvature index k,the 
gravitational held equations in GBD theory are 

d 2 + k a(f) toe /) 2 p . 

a 2 a(f) 6 (j) 2 3 (j) 


a dUfc d0 0 _ _P 

a a 2 2 <j) 2 a<fi 4> <f> 

where p and P are respectively the energy density and pressure of the fluid 
distribution.The equation of state of the fluid is given by P = yp.Some of the 
values of 7 for typical cases are -1 (vacuum),0(dust),l/3(radiation),l(massless 
scalar held). The wave equation for scalar held in the GBD theory is [2] 


4 > 



p — 3 P u ;0 

2u + 3 2u + 3' 


( 3 ) 


Energy conservation equation, which can be obtained from eqs.(l),(2) and 
(3) is, 


p + 3-(p + P)=0. (4) 

a 

To obtain a cosistent solution for the equations of motion we assume simple 
power law dependence on time for a ,(j) and uj in the following form which 
provides simple expanding solutions [ 6 - 10 ], 


a = ao 



( 5 ) 


* = (0 ■ (6) 

/ t \ H/3 

u((j)(t)) » (J) n = u 0 i^—J . (7) 

a is a positive number to ensure an expanding universe and n is a non zero 
constant. a,f3 and n are the free parametres at this level whose values we 
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wish to fix by way of obtaining consistent solutions of the gravitational and 
scalar field equations of motion. ao, 0o are present values of a(t) and 0(f).We 
further assume k = 0 to remain consistent with inflationary paradigm,which 
has received strong observational support.Before proceeding further,it may 
be noted that eqn(4) has the solution, 

p = Ca _3(1+7) 


and so 


P 


Po 


/ t \ -3+1+7) 

\TJ 


( 8 ) 


Po = Cciois the present(f = to) value of energy density. 

We now discuss the solutions of the field equations in different epochs. 


2.1 Radiation dominated epoch: 7 = 1/3 

For radiation dominated epoch the gravitational field equation [eqn.(l)]and 
scalar field wave equation [eqn.(3)] for k=0 become 


P 


a0 u) (0 
1 30 a0 6 l 0 


a\ 

aj 


and 


L 10 


+ 3-0 + 

a 2u + 3 


= 0 . 


Using equation (5),(6),(7) and (8), eqn.(9) yields 


(«Y = _po_ (±\- 4a -P __ M , (±\ np P 
V t J 3</>o \to J t 2 6 \io/ t 2 


n/3 


This equation can be put in the form, 


Pop 


(a + ap) - 3rPo 
This equation is satisfied when 


(J 


—4«—P+2 _ wo/3 2 


(+ 


n/3 


(9) 


( 10 ) 


—4a — ft + 2 = n(3 
a 2 + a/3 = 0 
02 _ _^Po/o 

Ct) O 0O 
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( 11 ) 

( 12 ) 

(13) 



Eqn.(lO), in similar manner leads to 


0 = o ( 14 ) 

0= (1 ~ 3tt )( 2 ^ + 3 ) ( 15 ) 

(n + 2)u 0 + 3 

Here we have further assumed 

fi(f) = 2c v(t) + 3 = f2 0 (^-) 13 , where = 2ca 0 + 3. 

Power law solutions (5),(6),(7) satisfy eqns.(9)and (10) simultaneously when 
eqn.(ll)and eqn.(12) are simultaneously satisfied with (3 given by either 
eqn,(14) or eqn.(15). 

The [3 — 0 solution of eqn.(14) implies a = 1/2 (from eqn.(ll)) and 
</>=constant and a;=constant which represents the well known Standard Model 
solutions.Thus, (3 = 0 solution represents the GR sector of GBD theory. 

Now using eqn.(15) in eqn.(ll)one gets 

_ 3n—2 cjq—3 

2no;o—2a;o+9n—3 

Similarly,eqn.(15)in conjunction with eqn.(12) gives 


2 ^o + 3 ( ^ 

a = --r- (lo) 

(4 - n) u 0 + 6 

For eqns(ll),(12) and (15) to be simultaneously satisfied , above two values 
of a should be equal. Imposition of this condition leads to 

3cuon 2 (2cUg T 9 cuo T 9)u T ( 6 ayj T 15aio T 9) — 0 

which admits two solutions like, n = — ^ 0+i '> and — ( 2cJ ° +3 ) . 

Substitution of these two values of n in eqn.(17),leads to a = 1/3 and 
a = 77+6 res P ec tively. Once again we find from eqn.(15)for a = 1/3, (3 
is equal to zero implying constant (j> and constant a; .We neglect this solution 
as it does not capture the characteristics of BD theory.Now for the second 
solution a = we find (3 = — 3/(u>o + 6 ) and n/3 = • Hence, the time 

dependences are obtained as 


Ci{t) — Oq 


t \ uq+6 

To) 


(17) 
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m = <h 7 )~" 0+ ' 

2cjq -(-3 

w (*) = (£) “ 0+6 

As can be seen from above ,the time dependence of a(t) <f>(t )and u{t) 
are completely fixed once lvq value is known, tuo needs to be fixed from other 
physical considerations. One can recover the standard model result for a(t) 
from above equations when cno = 0.Further it is easily seen that for ensuring 
deccelerated expansion in the radiation epoch(0 < a < 1) uq should be 
greater than —3.The values for (3 and n/3 are obtained to lie in the range 
—1 < (3 < 0 and — 1 < n/3 < 2 respectively. 


2.2 Matter dominated epoch: 7 = 0 

Using eqn.(5),(6) and(7), in eqn(l) we get 


(?r= vsr urOo (u 


—3a 


which can be written in the form 

2-3a-(3 


potp ( t_ 
30o \ 


H- -|- (x/3^ — 



(±\ nf> 

6 \t 0 J 


This equation is satisfied when 


2 — 3 a — /3 = n/3 

(18) 

c? + a/3 — 0 

(19) 

n2 _ 9 PQto 

</Wo 

(20) 


Similarly,using eqn.(5),(6),(7) in eqn.(3) one gets 


^/3(/3 — 1) + 3 oi(3 + 


n(3 2 UQ 

Qq 



—3a—n/3 


This equation is satisfied when 


(3 — 2 = —3 a — n/3 


( 21 ) 
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and 


P(P — 1) + 3 a (3 + 


n(3 2 uo\ 0 O 


Or 


+2 

hfl 


Po 

O 0 


With the help of eqn.(20),eqn.(22) reduces to 

( X + 2ojo+3 + 2(2^o+3) ) ^ + ( 3Q! “ l )@ = 0 - 

The above equation leads to,0 = 0 
and 

(l-3a)(4fj 0 + 6) 
(5 + 2n)u>o + 6 


( 22 ) 


(23) 


Inserting f3 = 0,eqn.(21) leads to a = 2/3.This is the standard model result 
representing GR sector of GBD theory. Using eqn.(23), eqn.(19) leads to 

_ 4^o+6 

(7-2n)u 0 +12 

Also for eqn.(23),eqn.(18) leads to 

_ (n+l)(4cjQ+6)-2[5ii>o+2iJon+6] 

3[(n+l)(4^ 0 +6) —(5o)o+2w 0 n+6)] 

For consistent solution of field equations,these two values of a should be 
equal.And we find 

2u 0 n 2 + [2cUq + 9cn 0 + 6 ] + [5u;q + 16cn 0 + 12] =0. 

Which admits two solutions like,n = — A+hAl arid ,n = — (u>o + 2). 

The first value of n is neglected as it once again gives ot — \ and (3 = 0 as in 
radiation dominated case. We ignore this solutions since it does not capture 
the time-varying characteristic of 0(f) and u(t) specific to the GBD theory. 
The second value of n leads to a = -A—- 

ojo+4 

Hence,the relevant solutions in the matter-dominated epoch are, 

a(t)=a 0 (^)^ (24) 

to 

^) = ^(£) “ 0+4 
2 (^p+ 2 ) 

(£) “ 0+4 

Once again it is found that the time-dependences of a(t),0(t), and u(t) are 
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entirely controlled by a>o,the present value of ca as was also the case in the 
radiation-dominated epoch.For deccelerated cosmic expansion uo need be 
greater than —2. That way even very large values of lvq required for agree¬ 
ment with solar system measurements[2] can also accomodated.The presently 
observed accelerated expansion too can be accomodated if —4 < ca 0 < 
—2.The values of (3 and n/3 are found to lie in the range — 1 < (3 < 0 and 
0 < n(3 < 2 respectively. 


2.3 Vacuum dominated epochs = — 1 


Using eqn.(5),(6) and (7),eqn.(l) reduces to 

( \ — rtn / + \ 4 rvl-t 


(*Y = _£o_ (t_Y P - 4. ^ (±\ nf> UY 

v t) 300 \to) t 2 6 \to J \tj 

This can be written as „ 

This equation is satisfied when 


0 2 UO f t\ 

np _ po% (±\ 2 ~ 0 


6 \toJ 

30o V *o / 


2-/3 = 

n/3 

(25) 

a 2 + a/3 

= 0 

(26) 

(3Y, _ 
6 

Pot 2 0 

30o 

(27) 


Similarly,using eqn.(5),(6),(7),eqn.(3) reduces to 

(W - 1) + 3«/3 + =£*) | (Y 2 = %(£)' 

This equation is satisfied when 

2 — f3 = n/3 

W-l)+3 ^+^ = <^§ 5 ; 

The above equation with the help eqn.(27) reduces to 




( 29 ) 



The (3 = 0 solution of this equation which leads to inconsistency with respect 
to eqns(25) and (28),we believe,is indicative of the impossibility of obtain¬ 
ing exponential standard model expansion within the power law assump¬ 
tions considered in this work. Continuing further,equation(28) and (29) give 

_ 3(n-l)-6aj 0 

U 3(2oj 0 +3)(n+l) 

However,equation(25) and (26)give 


a = 


ra + 1 


(30) 


Now these two values of a are equated for consistent solutions of equation(l) 
and (3).Hence,we get 


n = — (2ca 0 + 5) (31) 

This value of n leads to a = —and (3 = -Vtt 

UJ 0 +2 r' UIq+2 

Hence, the exact time dependence of a(t ),4>(t) and u(t) in case of vacuum 
energy dominated epoch are obtained as, 

a(() = «o(|-)"" +S (32) 

<P(t) = <k 

2ojq -(-5 

u ( f ) = (i) “ 0+2 

Once again, as in earlier two cases discussed, the time dependence is con¬ 
trolled by cjq, the present value of w.For accelerated cosmic expansion we 
must have —2 < wo < —1. In fact,closer the value of u to —2 faster is the 
accleration which can mimic inflation.Such a power-law inflation was indeed 
obtained by Mathiazhagan and Zohri,and La and Steinhardt [3] in their pio¬ 
neering works.Further in this epoch we obtain (3 and n(3 values to lie in the 
range — 00 < (3 < — 1 and 3 < n(i < 00 respectively for above stated range 
of Uq. 


3 Present value of 

In order to look for answer to the recently observed accelerated expansion of 
the universe,structure formation,coincidence problem etc.some authers have 
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studied cosmological models in Brans-Dicke theory or some modified version 
of this theory[10] .A small negative value of uj at t —> oo lying in the range 
—2 < cuo < —3/2 has been a recurrent result in these studies ,in contrdiction 
with large values of the order of 600 [2] obtained from solar system mea¬ 
surements. A possible solution of this contradiction,we believe, may lie in 
considering a non-constant coupling function cn(</>) as in GBD theory.Thus, 
the value of such a function can change with the cosmic time and,in the limit 
t —► oo, it could agree with local measured values. 

Further,as evident from our solutions,the time dependence of a(t),<j>(t) 
and u((f>(t)) in various epochs are entirly controlled by cno,the present value 
of u and not by the values of u> during the epochs considered. This brings out 
the importance of the value of ujq. For our results obtained in this work to 
be consistent with the results obtained by earlier authers[10] ,cno can safely 
lie in the range —2 < cn 0 < — 1 excluding the cn 0 = —3/2 value for which 
BD theory breaks down.Thus,we arrive at a comprehensive picture of cosmic 
evolution in the important epochs with the possibility of obtaining sufficient 
inflation in vacuum dominated era as well as accelerated expansion in the 
present epoch.Given the simple assumptions considered by us, we consider 
this a very interesting result which signifies the power of a GBD theory. 


4 VARIATION OF NEWTON’S GRAVITA¬ 
TIONAL CONSTANT G: 


As is well known, for theories with constant a;,the possibility of variations of G 
is very small.Consideration of arbitrary coupling function uj (0), however,opens 
the possibility of variations of G.ln the weak field limit Nordtvedt[2] found 
an expression for the observed value of the gravitational constant 


G{t) 


The variation rate of G is [2] 


/4 + M0)\ 

\3 + 2u(<f>)) 



G .</> 2c i> 

G = > + (3 + 2cu)(4 + 2u/) 


(33) 


(34) 
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where over dots represent time derivative. 

The present variation rate,therefore,is given by 




(to 


(3 + 2cno)(4 + 2loq ) 


( 35 ) 


b>r typical values of cuolike —1.9 or —l,we obtain from the above expression, 


~ 1.3/fo ~ 2 x 10 10 yr 1 . Thus,the variation of G is safely below the 
observational limit,which is 4 x 10~ lc Vr -1 [10,ll] obtained within the present 
formalism. 


5 CONCLUSIONS 

In this work,we have studied a homogeneous and isotropic cosmological model 
in the context of a generalised scalar-tensor theory. We have solved con¬ 
sistently both the gravitational held equation and scalar held wave equa¬ 
tion for matter,radiation and vacuum dominated epochs of the cosmic evo¬ 
lution by assuming power law time-dependence for the scale factor,scalar 
held and coupling parameter as extension of our earlier work[ 8 ].The solutions 
are generic which deserves attention.For all epochs,the solutions depend on 
one parameter i.e,on the present value of to. Interestingly,general considera¬ 
tions about cosmic evolution can constraint this value to lie in the range 
—2 < u 0 < —l.This result is quite consistent with other results[10] obtained 
from various other considerations.Also,we have shown that, when uq is chosen 
very close to -2 one can produce enough inflation according to wish,which can 
solve inflationary problems very easily in the standard model.While obtain¬ 
ing usual deceleration when the universe is radiation or matter dominated 
we also exhibit the possibility of accelerated expansion in matter dominated 
era for a slightly different value of ujq. 

We also find that the B-D scaler held <p{t) is always a decreasing function 
of time in all epochs with the decrease becoming very significant in vacuum 
dominated case.The time dependence of u{t) ,on the contrary,does not exibit 
such monotonic temporal behaviour. 

We have also examined the time variation of G(t)in this solution and 
found that the solution gives a safe time variation for the chosen range of Uq 
which gives a strong support to the results presented here. 
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